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1. Introduction 

In Ref. [|I| we applied the subtractive method of |0, Q to derive the NLO gluon-quark- 
induced contribution for the Drell-Yan process in a form suitable for use in a Monte- 
Carlo event generator. In this paper we extend these results to provide the angular 
distribution for the leptons, and, equivalently, the helicity density matrix for W and Z 
boson production. 

As described in [0, 0, ^, an event generator using this subtraction method generates 
two classes of events: One class is obtained from the leading order (LO) subprocess 
by showering the initial and final state quarks; the other class is generated by starting 
with a subtracted next leading order (NLO) subprocess. As in [|1|, we only treat the 
gluon-quark induced subprocesses, qg Vq' — > q'W (V = •y*, or This NLO 
subprocess is not necessarily suppressed by a factor of as if the gluon distribution 
function is larger than the quark distribution function. The quark-antiquark induced 
NLO corrections are complicated by the presence of soft divergences [Q, so we defer 
their detailed treatment in our method to the future. 
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The polarization state of the gauge boson (7*, or Z^) determines the angular 
distribution of the decay leptons in the lepton pair rest frame. The general structure 
of the angular distribution is given by nine helicity cross sections corresponding to 
the nine helicity density matrix elements for the gauge boson In the LO parton- 
level subprocesses {qq V ^ W) only the transverse polarization of gauge boson 
contributes to the helicity density matrix. In NLO, in general, all three polarizations 
of the gauge boson contribute, but the subtraction term, associated with a coUinear 
divergence, only involves the same transverse polarization as in LO. For the 7*, only 
the 4 parity-conserving terms in the cross section are nonzero @; but for and 
, with their characteristic V — A and V + A couplings, there are two other nonzero 
hehcity cross sections f^. 

The organization of this paper is the following. In Sec. 0, we describe the treat- 
ment of lepton pair production in an event generator. Then we compute the effect of 
combining the LO parton-level cross section with the order part of the shower; this 
will be needed as the subtraction term in the NLO hard cross section calculation. In 
Sec. ^, we carry out the subtraction from NLO helicity cross sections to get the angular 
distribution in the NLO differential cross sections. For and Z^ production, we also 
give the helicity density matrix. Finally, we summarize and discuss our results in Sec. 
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2. Monte- Carlo algorithms 

2.1 QCD improved parton model and inclusive lepton pair production 

Our calculations are based on the QCD improved parton model, which shows that a 
hard scattering process initiated by two hadrons is the result of an interaction between 
their constituents, namely quarks and gluons |^. 

We consider the process A + B^V + X—>-ll' + X in a. collision of two hadrons 
A and i? at a center-of-mass energy-squared s = {Pa + PbY- When the momentum 
fractions carried by the partons are not too small (roughly > 0.01), LO subprocesses 
{qq V W) are dominant, with the NLO subprocesses being suppress by a factor 
of as{Q'^)- But at the higher energies of the Fermilab Tevatron {^/s = 1.8 TeV), the 
momentum fractions are often smaller, and the gluon density can be larger than the 
quark and antiquark densities, so that gluon-induced NLO subprocesses like qg —>■ 
Vq' q'lV are not necessarily suppressed by as relative to the LO subprocess. 

The LO cross section da^'^iy) / dQ^dydVL for producing a lepton pair with rapidity 
invariant mass-squared and relative solid angle Vt = {6, (p) in the Collins-Soper 
frame |^ is obtained by weighting the parton level cross section with the parton distri- 
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bution functions (pdf's) fi/A{xa,Q^) and fi/sixh, Q"^] 



dQ'^dydn 



+ a{v)v^'^{e)f^>/A{xa, Q^)U/b{x,, g^), (2.1) 
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where 



Vqq>{9) = {1 + cos'^e + 2AqAicose), 

Vg,q{e) =Vq^in~6). (2.2) 

The parton-level cross section (y{V) and the coefficients Ag and Ai depend on the boson 
being produced: 

• v = Y- 

aiin = — X — A, = A = 0, (2.3) 
Ibvr dQ'^s 

where Cq is the quark charge in units of the size of the electron charge e. The 
summation is over (g, g') in {{u,u), {d,d), {s,s), (c, c)}, where we restrict our 
attention to quarks that can be approximated as massless in the production of 
W and Z bosons. 



V = W^: 



3 47ra2|y,|2 g2 
a{W ) = —— X — — ' X 



IGtt 9s 16< (g2 _ m2,)2 + M2,r2^, 



^ X ^v^G^lV^^f ^5(q2 _ M^)5(W^ ^ Iv) (2.4) 



167r 3 ^ 

■Aq = M = 1, : = -A = 1, 

where Vqqi is an CKM matrix element, s„, = sin^^, Gi;' is the Fermi coupling 
constant, and B{yV —>■ Iv) is the branching ratio T{W Iv)/T\y. In the second 
line of this equation, we have used the narrow width approximation for the decay 
of W boson. For the , the summation is over the cases (g, q') in {(m, c?), (c, s)}, 
while for the W~ we have (g, g') in {(c?, m), (s,c)}. 



aiZ"") = — X . X 



167r 9s 164c4 (Q2 _ M|)2 + M|r^ 



X 



-x/2Gp[g^ + gi]^5(Q2 - M|)5(Z° ^ /I) (2.5) 



IGtt 3 s 
_ 2gyg^ 2/yZ^ 
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where qa = T^, qv = — 2eqS^, = cos6'^, = sinOw, B{Z'^ 11) is the 
branching ratio for Z^, and the summation is the same as in the case of 7*. In 
the second hne, we have again used the narrow-width approximation. 

The fractional momenta assigned to the incoming quark and antiquark are equal to the 
following kinematic variables: 

(2.6) 

With full perturbative QCD corrections taken into account, Eq. ( p.l|) is replaced 
by the factorization formula, 

= S/''^''?^- (2,7) 

where the ^'s are the momentum fractions of the incoming partons, and ddij/ dQ'^dydVt 
is a suitably constructed parton-level hard scattering cross section, which depends on 
Xa/^i, Xb/^j, the vector boson mass Q, the solid angle fl, and on the renormaliza- 
tion/factorization scale through fj,/Q and as(/x). Now the sum is over all pairs of 
parton flavors (quarks and gluons). The formal domain of validity of Eq. (|2.7|) is the 
asymptotic 'scaling' limit, analogous to the Bjorken limit in DIS, where s,Q'^ ^ 00 
with Xa and Xb fixed. 

2.2 Parton-shower algorithm 

The initial-state shower algorithm for vector boson production used in PYTHIA or 
RAPGAP is described in (see also H). We use the method described in |l| to 
define the splitting variable z, rather than the s method of 0], since the treatment of 
NLO corrections is simpler. 

The part of the algorithm used in an event generator that concerns us is the fol- 
lowing: 

1. Generate values of y and Q^, then generate a direction (i.e., Q) for the leptons 
according to the LO cross section for lepton pair production, Eq. (pTT|). This gives 
the initial values of the parton and lepton 4-momenta. (These initial values have 
zero transverse momentum for the partons and hence for the lepton pair. After 
showering, the parton 4-momenta will be modified.) The direction Q gives polar 
angles 6 and (p that we choose to interpret in the Collins-Soper frame [Q. 
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2. For branching on hadron A side: Generate a virtuality Qf for the incoming 
quark i, the first longitudinal momentum splitting variable Za, and an azimuthal 
angle 0' for this branching. The distributions arise from the Sudakov form factor 



Here, Q^ax is normally set equal to for and Z^, and for 7* it is often fixed 
to be — the hard scattering scale. The Sudakov form factor is the probability 
that the virtuality of quark i is less than Qf. The branching on hadron B side is 
done similarly with first splitting variable z^. 

3. Iterate the branching for all initial-state and final-state^ partons until no further 
branchings are possible. 

4. Compute the actual values of the parton and lepton 4-momenta. 
2.3 Subtraction term from first initial-state branching 



In Sec. H, we will calculate the hard 
scattering cross section for gq-mdnced 
process (Fig. |I]). The gluon, of momen- 
tum Pa, comes from hadron A and the 
(anti)-quark, of momentum pb, comes 




from hadron B. The NLO contribution pjg^^^ gluon-quark collision. Gluon 

we want to calculate is to be accurate comes from hadron A and (anti)-quark comes 



when the following conditions are satis- from hadron B. 
fied. 1) The incoming gluon and quark 

have virtualities and transverse momenta-squared that are small compared with Q^. 2) 
The intermediate quark, of momentum pi, has a virtuality of order Q^. 3) The lepton 
pair has a transverse momentum of order Q. To avoid double counting, it is necessary 
to subtract the corresponding contribution obtained from the showering algorithm ap- 
plied to the LO parton level cross section. It is this subtraction term that we calculate 
in this section. After this subtraction, the combination of the LO and NLO term will 
also be accurate when the transverse momentum of the lepton pair or the virtuality of 
the intermediate quark decreases. 

The subtraction term is obtained by multiplying the LO parton level cross section, 
from Eq. ( |2.1| ), by the first order term in the expansion of the Sudakov form factor 

^The final-state showering is organized similarly to the initial-state showering. However, we will 
not need it explicitly in this paper. 




(2.8) 
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Eq. ( p.8|) in powers of as{Q'^)- It is made differential in the momenta of the particles 
involved, and the gluon-to-quark splitting kernel is selected. To match with the defini- 
tion of the NLO hard-scattering, it should be calculated when the initial-state partons 
are on-shell and have zero transverse momentum, and the final-state quark is on-shell. 
That is, we subtract the first-order cross section in the showering approximation: 




+ case: antiquarks (from g) come from A, quarks from B. (2.9) 

Note that in the summation in this equation, we still sum over q, q'. The gluon from A 
splits into a qq pair, and in the first term the quark from the gluon interacts with the 
antiquark from B. In the second term, on the last line, we have the antiquark from the 
gluon interacting with the quark from B. 

The splitting kernel here is for gluon quark + . ^{Qt) 

antiquark: P{z) = Pg^qq{z) = |(1 - 2z + 2z^). Note 
that because we are doing a strict expansion in pow- 
ers of as(Q^), the scale argument of the pdf is Q^. 
The function Ci{Ql) is a cut-off function that gives 

the maximum value of Qf. Normally it is 6{Q^ — 

, ^ 1 . rRi TfT ; 1 • c ,^ Figure 2: Collins-Soper frame, 

yf), but as discussed m y , a ditterent choice of the ^ , ^ , , , 

, r Pa and Pb lie m the x-z plane 

cut-oft function could minimize the number of nega- r i ^ ■ ^ r , 

oi the lepton-pair rest irame and 

five weight events generated generated from the sub- . , . , d i d 
^ ^ ^ z-axis bisects Fa and —rs- 

tracted NLO cross section. The general solution to 

deal with a smooth cut-off function is given in |^ for a non-gauge theory. We plan to 
apply this method to QCD in the future. 

Now, working in Collins-Soper frame ^ (see Fig. Q), we reconstruct the 4- vectors 
for the momenta g, p'^, pa and pb of the lepton pair, the outgoing quark, the incoming 
gluon and the incoming (anti)-quark. To be consistent with [Q, they obey the following 
requirements: 

1. The z-axis is given by the Collins-Soper definition for the angles of the leptons; 
it bisects the angle between Pa and —Pb, where the momenta are for the incoming 
hadrons. This is shown in Fig. |[ (Bengtsson, Sjostrand and van Zijl made a 
different choice of axes in the overall CM frame.) 

2. The incoming partons, pa and pb have momentum fractions C,a and C,b relative to 
their parent hadrons, in the sense of light-front components. 

3. pI = -Ql 
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4:. pI = pI = pf = 0, and Pa and pf, have zero transverse momentum in the overall 
CM frame. 

5. = Q^, where Q is the invariant mass of the lepton pair. 

6. y — - \n — ^ — - = - In — , is the exact rapidity of the lepton pair. 

2 Pa- q 2 Xb 



7. The first splitting variable on the A side is defined by 



Zn. = 



(2.10) 



Notice that in this definition, the splitting variable Za is defined in terms of the kine- 
matic variable Xa — e^Qj \fs^ rather than in terms of the fractional momentum of the 
intermediate quark. 

We define the Mandelstam variables at the parton level as usual: 

S = (Pa+Pbf = Ca^bS 
i = {Pb - = ^aXa- 



.S ^ -Ql 



U 



{Pa -qf = ibXb- 



Xa^b + Xb^a 
X^ -^"^ 



(2.11) 
(2.12) 



-s. 



(2.13) 



Xa^b + Xb^a 

These satisfy the constraint that s-\-i-\-u — Q'^ — Q. Then the partons' 4-momenta in 
CoUins-Soper frame are given by: 



P'b 



^aEA{l, sin 77, 0, cos 77), 
^bEs (1, sin?7, 0, - cos 77), 
^^^ = ^(1, 0, 0, 0), 



(2.14) 
(2.15) 
(2.16) 



p'r 



^p'^ + Pb-q' 



{[LEa + ibEB - Ql [iaEA + ibEs] sin?7, 0, [iaEA - ibEs] cos 77) , (2.17) 



Pv- 



— (1, sin 6* cos 0, sin 6* sin 0, cos 6*), 

— (1, — sin^cos0, — sin^sin0, — cos^), 
2 



with 



Ea = ^^/lT^e-y, Eb = ^VT+^e+y, 



sm77 



— r 



COST/ 



Qt 



ut 



(2.18) 
(2.19) 

(2.20) 
(2.21) 

(2.22) 
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Convenient variables for the above subtraction term are Q^, y, Q, C,a and Qi- 
However, they are not so convenient for the NLO matrix-element calculations. So we 
transform the cross section in Eq. ( p.9|) in terms of more convenient variables for a hard 
gluon-quark scattering: Q"^, y, Q , and ^h- Then the cross section after integrating 
over the branching angle 0' is: 

'^'''''-''^^^ E aiV) V^Ad)Ci{Q\)P{za) (2.23) 



dQ"^ dy dil d^a d^b 27r 



+ case: antiquarks (from g) come from A, quarks from B. 



3. NLO hard cross section 



3.1 Unsubtracted NLO term 



The angular distribution of the lepton pair produced in Drell-Yan process has been 
calculated in [|, |T0| for large Qt to order a^. The unsubtracted gluon-quark induced 



NLO scattering cross section corresponding to Fig. |1] is: 

?E HV)f-A^b,Q')ma,Q')H 

12 //OS I s\ ( m=4 



'^'^unsubtracted 1 (^) (^s{Q ) ^ ^nr\ £ ft r\2\ £ It r\1\ ZjU -\-L 



dQ^ dy dQ d^a d^b 27r s ^ 



Q\Q^ + g) 

s'^iXa^b + XbQ 



:i + COs2 5) + J2^f!iV)Gn 



m=Q 



+ case: antiquarks (from g) come from A, quarks from B, 



2^ 



2tt 

q,q 



E HV)f-A^b,Q')ma,Q 



2\ 



Cr (^aj ^b) Tjc/t c 



X 



where 



Xb-ib 

m=4 ~) 

;i + cos^ ^) + E Atiy)GAQ. 0) (s.i) 

+ case: antiquarks (from g) come from A, quarks from i?, 



^ + (3 2) 

—st ' 



A^''(V^) = 4Af (V) = Af{y) = 4Af{V) 
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iu[{Q^ + sy + iQ^-if 



{Q2-t){Q^-u)[(Q2 



+ {Q'-m 



(V^) = Af{V) 
Af{V) 



-Af = A,Ai 



Q^ + u^ + 2su 



ut 



-Af{V) = A,Ai 
A,Ai{l + A''i) 



\ 8(Q2 _ t)(g2 _ ^) (g2 _ 5)2 + (g2 _ £)2^ 

+ u2 _ 2QH 



(3.3) 
(3.4) 
(3.5) 



\ {Q^ - U){Q^ - i) (g2 - 1)2 + (g2 _ £)2 



(3.6) 



y^^" = A- 



\ {Q'-i){Q' 



2i{u + i) 



Go( 
G'i( 



\ iQ' - i){Q' - u) (g2 - sf + (g2 - iy 

^(l-3cos2^^), 

2 sin^ 9 cos 20, 
2^2 sin 26* cos 0, 
4^2 sin ^ cos 0, 
2 cos6', 



with s, i and m being given by Eqs. (|2.11| - ^T3D . Also: 

H (4a, = 7:^727 Tv^ XaXb[XaXb + 4a4fe), 



G'=(ea, e^) 



iiebiXa^b + XbiaY 
XbjXqXb + ^a^b)[{XaXby + {XgXb - ^a^bf] 
^l^bi^a^b + Xb^a)iXb + 6) 



(3.7) 
(3.8) 

(3.9) 
(3.10) 
(3.11) 
(3.12) 



(3.13) 
(3.14) 



as defined in |TT| and used in [||]. 

Notice that, for V = 7*, Ag = Ai = 0, so (7*) = (7*) = 0, that is, only the 
first 4 helicity structure functions are nonzero for 7*. 

In the above Eqs. ( p.3| - p.6| ), only the term containing A'^'^{V) (V = W^, Z^) will 
give a collinear divergence. To get the hard cross section, we need to cancel out the 
divergence using the subtraction term obtained in Eq. ( |2.23| ), so we rewrite our NLO 
unsubtracted parton level cross section Eq. ( p.l[ ) as: 

'^0"unsubtractedl(^) C^siQ"^) 



dQ"^ dy dQ dC,a d^b 



2t\ 



E '^(^)/^'(6,g')/.(ea,g' 
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X 



X 



m=3 

+ case: antiquarks (from g) come from A, quarks from B. (3.15) 
3.2 NLO term with subtraction 

We now subtract the showering term, Eq. ( p.23| ). Only the term proportional to Vgg'{6) 
and T>-qiq{9) needs a subtraction term; all other terms are finite in the coUinear limit 
(t — 0). So we have: 



, (New) 
"'^hardl 



{V) 



dQ"^ dy dVt d^a d^b 



27r 



Xb - 6 



+ H'Ua,^b) 



f-M,Q')Ci{Ql)P{z,) 



'^Xbiaib + 



+ f-A^b,Q 



X 



Xb - 6 

E Af!{V)GUO, 0) + A,AiA'^G4 



m=3 



m=0 



2tt 



+ case: antiquarks (from g) come from A, quarks from B 

q,q' 



(3.16) 



+ 



G'iia, ib) 



X 



Xb - ib 

m=3 

E ^'i(y)G 

m=0 



+ H\ia.^b) 



Xb - ib 

+ case: antiquarks (from g) come from A, quarks from B 



G 



(3.17) 



where 



fjU _ G'^{^a, C,b) . , 
^ _c + ^ ^b) 

•^b So 
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Observe that, the coefficient of T>qqi{6) agrees with the results in detailed cal- 
culations shows that the pdf used here is the same as in []T| and [Q. 

3.3 Results for qg subprocess 

Now we specify the changes needed for the other gluon-quark scattering subprocess, 
in which gluon comes out of hadron B instead of A. The formulae for the subtracted 
cross section are obtained from Eq. ( p.l6|) by making the following substitutions in the 
right hand side of those equations: 

• Case that quark comes from A, antiquark from B: 

a ^5, 1^2, (V) ^ A^^(r). (3.19) 

• Case that antiquark comes from A, quark from B: 

a ^6, 1^2, A^^iV) ^ A?SV). (3.20) 

Here 

A^Ji = Af!{u ^ t), At = A^^iu ^ £), (m = 0, 1, 4); 

Al^ = -Af (^ ^ t), Al^ = -Ai\u ^ i), {n = 2, 3). (3.21) 

The cut-off function Ci{Ql) is replaced by a function C2{Ql) with the same functional 
form, and we define 

-Ql = u = sx,^, f ~ , z, = ^. (3.22) 

3.4 Helicity density matrix for and resonance 

In some cases, we want to separate the contribution of different spin configurations in 
the cross section, and then the helicity density matrix is very useful. Instead of giving 
the explicit formula for helicity density matrix, we give a prescription to obtain the 
helicity density matrix from the helicity cross sections. 

In Eq. ( |3.1(j| ), inside the curly brackets, let us define , H^, H^, H^, and 
to be the coefficients of the angular distributions (1 + cos^6'), Go{9,(j)), Gi{9,(j)), 
G2{0,(f)), G-i{9,(j)) and 6*4(6', 0). In [|], these are expressed in terms of what are called 
'helicity density matrix elements', but the normalization condition of a density matrix. 
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that the trace is unity, is not satisfied. So we will simply called them 'helicity matrix 
elements'. Thus the if's are linear combination of the helicity matrix elements. The 
helicity matrix H^l^' is then given by 



ttOO 
Hab 


= Hab 






^ab 


= {Hi 


- Hab 


+ H^b)l'^ 


Hab' 


= (Hab 


- Hab 


- ^r.)/2 


Hab~ 


= Hab^ 


= Hlb 




^ab 


- J^ab 


= Hab 




Hab 


= Hab 


= Hab 


^ab^ 



(3.23) 

where the subscript ab labels the parton content of the cross section. Then the actual 
helicity density matrix is obtained from the above matrix by dividing it by if^. 

In the subprocess gq Vq' —>■ q'll', we have ab = gq, while in the subprocess 
qg Vq' q'W, we have ab = qg. 

4. Conclusion 

We have extended the application of Collins's MC algorithm to give the angular distri- 
bution of Drell-Yan lepton pair production at order a^. We also gave the corresponding 
helicity density matrix of W and Z bosons. We confirmed that the pdf 's used in Collins 
algorithm are process- independent at order a^. 
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